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ABSTRACT 

We present A^-body simulations of dissolving star clusters close to galactic centres. 
For this purpose, we developed a new A^-body program called nbody6gc based on 
Aarseth's series of A^-body codes. We describe the algorithm in detail. We report about 
the density wave phenomenon in the tidal arms which has been recently explained by 
Kiipper et al. (2008). Standing waves develop in the tidal arms. The wave knots 
or clumps develop at the position, where the emerging tidal arm hits the potential 
wall of the effective potential and is reflected. The escaping stars move through the 
wave knots further into the tidal arms. We show the consistency of the positions of 
the wave knots with the theory in Just et al. (2009). We also demonstrate a simple 
method to study the properties of tidal arms. By solving many eigenvalue problems 
along the tidal arms, we construct numerically a ID coordinate system whose direction 
is always along a principal axis of the local tensor of inertia. Along this coordinate 
system, physical quantities can be evaluated. The half-mass or dissolution times of our 
models are almost independent of the particle number which indicates that two-body 
relaxation is not the dominant mechanism leading to the dissolution. This may be 
a typical situation for many young star clusters. We propose a classification scheme 
which sheds light on the dissolution mechanism. 
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1 INTRODUCTION 

The centres of the Milky Way and othergalaxies are cur- 
rently a field of very intensive researchrj In our Galaxy, 
observations have to be carried out in other wavelengths 
than visual due to the huge extinction. Directly in the cen- 
tre of our Galaxy resides the strong radio source Sgr A* 
at the location of the Galactic super-massive black hole 
(M, « (3 - 4) X 10*^ Mq, e.g. Genzel et al. 2000, Ghez 
et al. 2000, Schodel et al. 2002, Ghez et al. 2003, Eckart et 
al. 2005, Ghez et al. 2005, Beloborodov et al. 2006). The 
Galactic centre region spans roughly nine orders of mag- 
nitude in galactocentric radii ranging from a rough outer 
radius of the central molecular zone (i?cMZ ~ 200 pc, Mor- 
ris & Serabyn 1996) down to the Schwarzschild radius of the 
Galactic super-massive black hole (ii. 4 x 10"'^ pc). This 
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large range in radial scales already suggests that the physics 
in the Galactic centre region is extremely rich in content. 

Two young star-burst clusters named Quintuplet (Na- 
gata et al. 1990, Okuda et al. 1990) and Arches (Nagata et 
al. 1995) have been discovered at projected distances less 
than 35 pc away from the Galactic centre. They have quite 
extraordinary properties and stellar contents. Their forma- 
tion still requires clarification. However, both clusters are 
located (at least, in projection) near the Galactic centre 
"Radio Arc" (Yusef-Zadeh, Morris & Chance 1984, Tim- 
mermann et al. 1996), which is a region rich in molecular 
clouds and gaseous filaments (Morris & Serabyn 1996, Lang 
et al. 2005). 

The tidal field is extremely strong in the Galactic cen- 
tre. It was therefore highly desirable for us to study the eifect 
of the tidal field on the dynamics of star clusters which orbit 
around galactic centres at small galactocentric distances. A 
few similar simulations as those presented in this paper can 
be found in the works by Fujii et al. (2007, 2009). Other pre- 
vious works on the dissolution of star clusters in the Galactic 
centre have been published by Portegies Zwart, McMillan & 
Gerhard (2003), Kim & Morris (2003) and Guerkan & Rasio 
(2005). 
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Figure 1. Sketch of the geometry of the problem as seen from 
the galactocentric reference frame. The vector points from the 
Galactic centre to the star cluster centre. The vector points 
from the star cluster centre to the position of the ith star. The 
orbital velocity of the star cluster has been denoted as Vg. 



We do not attempt to solve the paradox of youth (Ghez 
et al. 2003) in this study with the star cluster in-spiral 
scenario (Gerhard 2001). Gerhard used In A = 10 for the 
Coulomb logarithm of dynamical friction. This value is, from 
our point of view, much too large. We will use in this study 
a more realistic and variable Coulomb logarithm according 
to Just & Pefiarrubia (2005) which leads to in-spiral time 
scales which are considerably larger. 

This paper is organised as follows: In Section 2, we 
describe in detail the algorithm of our A'^-body program 
nbodyGgc which has been especially developed to study 
the dynamics of star clusters in galactic centres. Section 3 
describes the theoretical models which we use for the cen- 
tral region of the Galactic bulge and the star clusters. In 
addition, we discuss the effective potential which is essen- 
tial in order to understand the dynamics in the tidal field, 
and Poincare surfaces of section. Section 4 contains the re- 
sults of our direct A'^-body simulations. Our main focus is 
on the properties of the tidal arms and the dissolution time. 
In Appendix [X] we show the Taylor expansion of the effec- 
tive potential. In Appendix |B] we describe the algorithm of 
an eigensolver which is used to construct a ID coordinate 
system along the tidal arms. It can be used to evaluate phys- 
ical quantities along the tidal arms in order to study their 
properties. 



2 NUMERICAL METHOD 

The computer program nbodyGgc which is used in this 
study is a variant of the A'^-body program nbody6+ + 
(Aarseth 1999, 2003, Spurzem 1999) suited for massively 
parallel computers]^ The code nbodyGH — h is a variant of 
the direct A^-body code nbody6 (Aarseth 1999, 2003) for 
single-processor machines. A fourth-order Hermite scheme, 
applied first by Makino & Aarseth (1992), is used for the 
direct integration of the 3A'^ Newtonian equations of mo- 
tion of the A'^-body system. It uses adaptive and individual 

We remark here that NB0DY6GC is based on a code variant 
called nbody6tid which has been developed by R. Spurzem 
in collaboration with O. Gerhard and K.-S. Oh (unpublished). 
nbody6tid was very helpful for the development of NB0DY6GC. 
However, we switched to another integrator for circular and very 
eccentric cluster orbits and improved the treatment of dynamical 
friction for studies in the Galactic centre. 



time steps, which are organised in hierarchical block time 
steps, the Ahmad-Cohen neighbour scheme (Ahmad & Co- 
hen 1973), Kustaanheimo-Stiefel (KS) regularisation of close 
encounters (Kustaanheimo & Stiefel 1965) and Chain regu- 
larisation (Mikkola & Aarseth 1990, 1993, 1996,1998). 

2.1 Cluster orbit 

We denote the radii, velocities and accelerations related to 
the Galactic centre with capital letters and those related to 
the star cluster with lower case letters. Figure [l] shows the 
geometry of the problem: A star cluster is orbiting around 
the Galactic centre. The potential in which the star cluster 
moves is the sum of the Kepler potential of a super-massive 
black hole and a scale free pote ntial of the central region 
of the Galactic bulge (cf. Section 3.11 ^The two first-order 
equations of motion for the star cluster orbit read 



V9(t) = Rj, 

V,{t) = -V<1>,(|R,1) + Ad 



(1) 
(2) 



where R9,Vg,$g and aaf are the position vector, velocity 
vector, gravitational potential of the Galactic centre region 
and deceleration due to dynamical friction and the dot de- 
notes the derivative with respect to time. The equations of 
motion ([l]) and ([2| for the star cluster orbit with respect to 
the Galactic centre are solved using an 8th-order composi- 
tion scheme (McLachlan 1995; for the idea see Yoshida 1990) 
with implicit midpoint method (e.g. Mikkola & Aarseth 
2002), thereby including a realistic dynamical friction force. 
Although the symplectic composition schemes are by con- 
struction suited for Hamiltonian systems, they can be used 
for dissipative systems as well if the dissipative force is not 
too large. In our case, four iterations turned out to be suffi- 
cient to guarantee an excellent accuracy of the scheme. 

We use a cluster membership criterion such that the 
dynamical friction force is based only on the total mass of 
the cluster members. We define a membership radius rm by 
the condition 



Pel = 



3Mci(r. 



(rm) 



(3) 



as the radius where the mean density p^,; in the star cluster is 
equal to the local bulge density at the cluster centre which 
is located at radius Rg. This radius differs from the tidal 
radius (King 1962) only by a factor of order unity. Stars 
within twice the membership radius are defined as cluster 
members. 



2.2 Stellar orbits 

On the other hand, the equations of motion for the or- 
bits of stars in the star cluster are solved by the standard 
nbody6/nbody6++ routines using the 4th-order Hermite 
scheme (Makino & Aarseth 1992), KS or chain regularisa- 
tion including the full 3D tidal forces from the super-massive 
black hole and the Galactic bulge. The tidal force is added 

^ The program nbody6gc is written in a way that any analytical 
galactic potential can be implemented. 
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as a perturbation to the KS regularisation. The following 
quantities are involved: 

(i) The specific force on the ith particle due to all other 
stars (cluster members and non-members) is given by 



(4) 



where A'^, G, rji = Vj — r^, rrij are the particle number, the 
gravitational constant, the relative position vector between 
the ith and jth particles and the mass of the jth particle, 
respectively. 

(ii) The specific force due to the Galactic centre at the 
position of the cluster centre is 



GC GMbh 



Rg " Rg J Rg 



R.IJ 



(5) 



where C, Mhh and a are the normalisation of the scale free 



bulge mass profile (see Section 3.1 1, the mass of the super- 
massive black hole and the cumulative mass profile power 
law index. 

(iii) The specific force exerted on particle i due to the 
Galactic centre is given by 



A.gi 



GC ^ GMbh \ Rgi 

R^T" Rni I Rfli 



(6) 



by 



(iv) The deceleration due to dynamical friction is given 



AnG'^PgMci 



\nKx{Vg) 



(7) 



where pg, Aid, Vg and Vg are the local bulge density at 
the position of the star cluster centre, the star cluster mass 
and the velocity vector and modulus of the Galactic centre, 
respectively. Furthermore, In A is the Coulomb logarithm 
which results from the integral over impact parameters and 
xi'>^g) = /o^' f{''^)d^''^ is the result of the integration of the 
distribution function f{v) of light particles over velocity 
space. For the Coulomb logarithm In A, we use according 
to Just & Peiiarrubia (2005) 



In A 



1,2 



bo 



rv, 



Pg 



(8) 
(9) 



where fei, bo, L are the maximum and minimum impact pa- 
rameters and the local scale length of the bulge density pro- 
file, respectively, and rv ~ GM^i/{i\Eci\) ~ is the virial 
radius of the star cluster (where Ed is the internal energy 
of the star cluster and is the half-mass radius). 

In the galactocentric reference frame, the total force on 
the ith particle would be given by 



*-tot,i,gc 



*-tot,i,gc 



— <^i.qc ~t- ^-qi qc 



memb. (10) 
non — memb. (11) 



where the subscript "gc" denotes "galactocentric" . However, 
we choose the cluster rest frame as reference frame for our 
simulations. This is necessary, because Aarseth's family of 
Ai'-body programs is adapted to this reference frame and 
assumes that the cluster centre is close to the origin of co- 
ordinates. This guarantees a sufficient accuracy of the Her- 
mite scheme which is used for the orbit integration. On the 
other hand, this choice of the reference frame implies that 
the Galactic centre is modelled as a pseudo-particle which 
orbits around the cluster centre. We keep in mind that a 
transformation from the galactocentric frame to the clus- 
ter rest frame implies that r^, Rg, Vg and Rgi in Q - ([7| 
change their sign. This implies that 



^i,cl *^i,gC5 ^g,cl ^-g.gcj (-i^) 

^-gi,cl — ^-gi,gc^ ^-df,cl — ^-df,gc (-i"^) 



where the subscript "cl" denotes the cluster frame. It is then 
convenient for the force computations to transform to a ref- 
erence frame in which the initial cluster centre is force-free. 
Since this frame is accelerated, an apparent force 



^g,cl ^df ,cl 
-■^g,cl 



memb. 

non — memb. 



(14) 
(15) 



appears. In the accelerated cluster frame the total force on 
the ith particle is therefore given by 



■^tot.i.acl — -A-to; 



t,i,gc 



(16) 



where the subscript "acl" denotes the accelerated cluster 
frame. Thus the second-order equations of motion for the 
orbits of the cluster stars read 

SLtot,i,aci = a-i.ci + Agi,ci — Ag,ci mcmb. (17) 

9.tot,i,acl = a^.c! + ■Agi,cl — ■Ag,cl — Adj\cl n. — m. (18) 

It can be seen that in the accelerated cluster frame an in- 
dividual star experiences only the differential tidal force be- 
tween its own location and the cluster centre. 

We applied a density centre correction in certain inter- 
vals to correct for the displacement of the density centre. 
This was done in order to retain a consistent treatment of 
dynamical friction since the dynamical friction force is de- 
termined from the approximation that the star cluster mass 
is concentrated in the origin of coordinates. 



2.3 Energy check 

The specific energy of a particle is calculated in the galac- 
tocentric reference frame: 



1 f'^ 
e,; = -(vi -I- Ygf $i,,„t + -l-^.e^t + / ad/ • Ygi dt (19) 

^ J to 

where $i,i„t = — GJ^^^^.^-j mj/{\vi — rj\) is the fuU inter- 
nal potential of the iV-body system, '^i,ext is the external 
potential of the Galactic centre and the last term is the en- 
ergy loss due to dynamical friction. Using these terms, we 
calculate a total energy E\. A factor 1/2 has to be included 
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in the summation of the potential energy of the A'^-body sys- 



tem. Note that the individual terms in ( 19 1 have quite dif- 
ferent orders of magnitude. On the other hand, the orbital 
energy E2 of the star cluster in the galactocentric reference 
frame is given by 



E2 



a 



■f 

Jo 



Mcl ad/ 



Vgdt + 



* Mci 



Eodt 



(20) 



where the last two terms on the right-hand side are cor- 
rections due to dynamical friction and tidal mass loss of 
the star cluster, respectively. Both energies Ei and E2 are 
checked at regular intervals for conservation. We note that 
the problem of a star cluster orbiting in a galactic tidal field 
has two energy scales related to the internal energy of the 
star cluster and the external energy of the tidal field. Tidal 
heating can transfer external energy from the larger scale 
into internal energy. The program nbody6gc conserves the 
energy within the cluster to a sufficient degree such that 
two-body relaxation is not suppressed (cf. Ernst 2009 for 
more details). 



Parameter 


Value 


Parameter 


Value 


a 


1.2 


Ej{Li) [pc2Myr-2] 


1.67546e5 


Ro [pc] 


20 


Ej{L2) [pc2Myr-2] 


1.67592e5 


Mo [Mq] 


1.67459c8 


Ej(L3) [pc2Myr-2] 


1.65965e5 


Po [M0pc-»] 


1998.90 


<S>.a.tid{Ro) [pc2Myr-2] 


1.69502e5 




1.88335e5 


x(Li) [pc] 


-2.66618 


ojo [Myr-i] 


9.704 


x{L2) [pc] 


2.78522 


Ml [Mq] 


10® 






ri [pc] 


1.20213 


G [pc3 Mg^ Myr-2] 


(222.3)-^ 



Table 1. Parameters used for the model of the Galactic centre 
region (Section |3.1| and the Plummer models (Section |3.2[ l which 
are used in Sections |3.3| and |3.4| C{Li) and x{Li) the value of 
the effective potential at the Lagrange point Li and its location, 
respectively and G is the gravitational constant. 







dhiLO 
















din It 





+ 1 



(27) 



The angular momentum of the circular orbit is given by 



3 THEORY 

3.1 Galactic centre model 

In the following Sections, we will use parameter values close 
to those of the centre of the Milky Way, since these param- 
eters are better known than those of any other centre of a 
galaxy. 

For the very central region of the Galactic bulge, we 
use a spherically symmetric scale free model (i.e. a model 
which is self-similar under scaling of lengths). The potential 
"1>, cumulative mass M and density p are given by 



M{u) 
p{u) 

where 



<I>o In u a = l 

Mou" 

3 



R/ Rq, 



(21) 

(22) 
(23) 



Po 



a Mo 



$0 



vPoRo, 



(24) 
(25) 



1 GMo _ AnG 
a — 1 Ro a{a — 1) 

a is the power law exponent of the cumulative mass pro- 
file, G is the gravitational constant and Ro is a length unit 
(which is not inherent in nature but simply a human con- 
vention) . 

The circular frequency uj is given by 



/ \ {q-3)/2 / 



47rGpo 



(26) 



The ratio of the epicyclic frequency k to the circular fre- 
quency Lu is given by 



L{u) — Lou 



(c+l)/2 



Lo — uJoRo- 



(28) 



The parameters of our models are given in Table [T] The 
value of Mo corresponds to Mo{Ro = 1 pc) = 4.6 x 10^ Mq. 
The stellar mass within the central parsec is not easy to 
determine (see Schodel et al. 2007, Genzel et al. 2003 and 
also the review by Mezger, Duschl & Zylka 1996). 

The simple numerical calculations in Sections |3.3| and 
13.41 have been done without a black hole at the Galactic 
centre. Nevertheless, in our A'^-body calculations, we added 
the contribution of a super-massive black hole of mass M, — 
3.6 X 10® Mq (Eisenhauer et al. 2005). However, the influence 
radius (Frank & Rees 1976) of the Galactic super-massive 
black hole is only 1 — 2 pc which is small compared to the 
galactocentric radii used in this study. 

3.2 Star cluster model 

For the star clusters, we use Plummer models for the sim- 
ple numerical calculations in Sections |3.3| and |3.4| and King 
models (King 1966) for all A'^-body models in Section [4] 
The Plummer model is given by 



$(w) — — <I>i 
M{v) = Ml 
p{v) 



[l-h«2]3/2 
1 



1 + v 



215/2 



(29) 
(30) 
(31) 



with the dimensionless radius v = r/ri ^ and the Plum- 
mer radius 



ri 



_ GMi _ / 3Mi \ 
$1 \4:npi J 



1/3 



(32) 
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Figure 2. Effective potential for a star cluster in the centre of a galaxy {z = plane). The large potential well is due to the Galactic 
centre and the small one is due to the star cluster. 



King models {Wq = 6): 

In general: A/ci = W^Mq^t^ = 1.64 pc, n = 11.2 pc; 
Kl (N = 103), K2 (Af = 2 X 10^), K3 (A^ = 5 x lO^), 
K4 (N = 10*), K5 (Af = 2 X 10*), K6 (Af = 3 x 10''), 
K7 (Af = 5 X lO-*), K8 (Af = 7 X lO-*), K9 (A^ = lO^) 



Table 2. Parameters of the TV-body runs with King models. Wo, c 
and rfi are the dimensionless central potential, the concentration 
and the half-mass radius of the King model, respectively, A'^ is the 
particle number and M^i is the total cluster mass (which is only 
needed to calculate a dynamical friction force). 

where Mi is the total cluster mass, — $i is the central po- 
tential and pi the central density, all of them being finite. 
The parameters of the Plummer models are given in Table 
[l] The parameters of the King models which we use for the 
A'^-body simulations are given in Table [2] 




3.3 Effective potential 

In this and the following Section, we use the parameters 
given in Table [l] The length unit _Ro corresponds to the 
radius Rc of the circular orbit, i.e. we have Ro = Rc- 

Since we are considering a circular orbit, it is convenient 
to study the physics in a reference frame which is co-rotating 
with the frequency luc = of the circular orbit. The star 
cluster centre is taken as the origin of coordinates. We choose 
a right-handed coordinate system where the a;-axis points 
away from the Galactic centre and the y-axis points in the 
orbital direction of the star cluster orbit around the Galactic 
centre. In this reference frame, centrifugal and Coriolis forces 
naturally appear according to classical mechanics. The po- 
tential in which a particle moves is the superposition of the 



effective tidal potential and the star cluster potential. For 
short we will call this the effective potential. 

The effective potential is shown in Figure [2] It is given 
by the expression 



^ix + Ro)^+y2 + z'^ 
Ro 



^ujo [{x + Ro f + y'^ + z^] 
GMi 

sjRl + x'2 + y2 + z2 



(33) 



Note that $o, <^o and _Ro are related by <l>o = LJoRo/{a — 1). 



In Equation ( 33 1, the first term is the gravitational potential 
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tidal potential $off,tid- The corresponding ID Taylor series 
of "I'cff ,tid along the x-axis around x — {y = z — 0) is given 
by the power series 



1.69 



1.68 - 



•B. 1.67- 




1.66 = 1 



Figure 4. Zoom into the star cluster region of Figure [2] along 
the X axis with y = 0. The Galactic centre lies in negative x 
direction. It can be seen that the Lagrangian points Li and L2 
lie at different energies and at different distances from the cluster 
centre due to the asymmetry of the effective potential with respect 
to a: = 0. The dashed line marks the effective tidal potential. The 
solid line is the full potential with the added contribution of a 
Plummer potential. 



of the central region of the Galactic bulge, the second term 
is the centrifugal potential and the last term is the Plummer 
potential of a star cluster. Note that the bulge potential is 
not well behaved in the limit Ro ^ if there is no black hole. 
However, the physics considered in this work happens close 
to the radii of the circular orbits. We note that the Jacobi 
energy per unit mass Ej — (i^ + + i^) /2 + "l>cff(2;, y, z) 
is a conserved quantity in the co-rotating reference frame. 
The tidal terms (i.e. the first two terms on the right- 



hand side) of Equation ( 33 1 can be expanded in a Taylor 
series around the star cluster centre {x,y,z) — (0,0,0). Up 
to the 5th order, the solution is given in Appendix [A] The 
expansion up to the second order coincides with the tidal ap- 
proximation which is a linear approximation of tidal forces. 
This approximation can be used to study the dynamics in 
star clusters on circular orbits which are far away from the 
Galactic centre. We stress, however, that we used the exact 
expressions for all computations in this study. 

Figure |3] shows a zoom into the equipotential lines 
around the star cluster region of Figure [2] Figure |3] also 
shows the location of the Lagrange points Li and L2. As 
usual, Li lies on the negative i-axis (between the cluster 
centre and the Galactic centre) while L2 lies on the positive 
a;-ax;is. Li and L2 are saddle points of the effective potential. 
It can be seen that at the locations of Li and L2, the surface 
in Figure [2 is curved differently along the x— and j/— axes. 

Figure 4]shows the effective potential in the star cluster 
region along the line connecting the Galactic centre with the 
star cluster centre. The dashed line shows only the effective 



$off,tid 



ujqRo + o(a - S)uJoX 

a ~ i / Z 



n(a-o 



2 



(34) 



Higher-order terms lead to an asymmetry with respect to 
X — which becomes important in the vicinity of the Galac- 
tic centre. A non-linearity in the tidal forces is related to this 
asymmetry. Such non-linear effects can be seen in Poincare 
surfaces of section. The solid line in Figure |4] shows the full 
effective potential. It is the superposition of the effective 
tidal potential and the star cluster potential. The Lagrange 
points Li and L2 lie at slightly different energies and at 
slightly different distances from the star cluster centre whose 
position we denoted as as L3. The energies and locations of 
the Lagrange points are given in Table [l] 

We stress that this picture is only valid for a star cluster 
orbit which is exactly circular. The region above the tidal 
effective potential (dashed line in Figure [4| is energetically 
forbidden for the cluster orbit. As soon as it becomes eccen- 
tric, the cluster centre no longer remains at the position of 
the extremum of the effective tidal potential but oscillates 
around a:: = and is reflected either at the centrifugal or the 
gravitational barrier. This oscillation leads to oscillations of 
the Jacobi energies of the Lagrangian points Li and L2 on 
the orbital time scale of the star cluster orbit and can change 
the dynamics dramatically. 



3.4 Poincare surfaces of section 

Figure [5] shows a few Poincare surfaces of section for the 
orbit with Parameters given in Table [T] which is exactly cir- 
cular. 

The left column of Figure [5] shows two Poincare sur- 
faces of section at a Jacobi energy deep in the potential well 
of the star cluster. The equipotential line corresponding to 
this Jacobi energy (which corresponds to the envelope of the 
lower surface of section in the left column of Figure [5f al- 
most has a circular shape. The Poincare surfaces of section 
at this Jacobi energy show that all orbits are regular and 
confined to invariant curves by a third integral. Such a third 
integral can usually be represented by a power series expan- 
sion where the lowest order is the angular momentum which 
would be exactly conserved if the system were spherical. 
Since the system is in fact not exactly spherical, the angu- 
lar momentum slightly oscillates around some value (see e.g. 
Figure 3-5 in Binney & Tremaine 1987). 

The middle column of Figure [5] shows two Poincare sur- 
faces of section at the Jacobi energy Ej = Ej{L2) which 
corresponds to the Lagrange point L2. The equipotential 
lines are open around Li and particles can escape towards 
the Galactic centre. The phase space is divided between reg- 
ular and chaotic regions. 

The right column of Figure [5] shows two Poincare sur- 
faces of section at a Jacobi energy which is higher than the 
value of the effective potential at both Lagrange points Li 
and 1/2- The equipotential lines are wide open around Li and 
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Figure 5. Poincare surfaces of section. Left column: Deep in the potential well of the star cluster at C = 1.66638e5 pc^/Myr^. Middle 
column: At Ej = Ej{L2)- Right column: Above the Jacobi energies of Li and L2 at Ej = 1.68845e5 pc^/Myr^. 



L2 and particles can escape in both directions either into the 
leading or the trailing tidal arm. All orbits are chaotic. 



4 RESULTS 

4.1 Tidal arm properties 

Figures [6] -[9] show the formation of the tidal arms for model 
K9. The initial 90% Lagrangian radius has been taken to 
be equal to the membership radius in Equation Q. The 
star cluster dissolves in a spiral-like structure. The leading 
tidal arm consists of particles which pass the inner Lagrange 
point Li , while the trailing arm is formed by particles which 
pass the outer Lagrange point L2 . The galactocentric radius 
Rc of the circular orbit is shown as a solid line. It decays 
very slowly due to dynamical friction with the Coulomb log- 
arithm which was modified according to Equations ^ and 
([9|. The initial value for the model K9 is InA 1.7. Most 
particles of the leading arm have galactocentric radii less 
than Rc while most particles of the trailing arm have radii 
larger than Rc- The dashed lines mark once and twice the 
membership radius r™. 

We have introduced a local coordinate system according 
to the description in AppendixjB] We denote the coordinate 
along the tidal arms as w, where negative values refer to 
the leading arm and positive values to the trailing arm. The 
short and long marks correspond to multiples of 1 and 5 pc, 
respectively. 

The color coding is according to the logarithm of the 



0,50 



log p [Ma/pc' 



6,50 




-40 L 



-+0 



-20 







20 



40 



Figure 6. The model K9 at t = 0.22 Myr. The Galactic centre is 
marked with a cross. The star cluster orbit is shown as a solid line. 
The dashed lines mark once and twice the membership radius. 
We look in the direction of the Galactic north pole. The short 
and long marks of the tidal arm coordinate system correspond to 
multiples of 1 and 5 pc, respectively. 
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Figure 7. Further evolution of the model K9. Top panel: At 
t = 0.43 Myr. Bottom panel: At t = 0.87 Myr. The dotted hues 
from the Galactic centre show the angles between density maxima 
in the leading and trailing arm, respectively, with respect to the 
star cluster centre. 

stellar density. The density clearly peaks in the cluster cen- 
tre. However, one can observe clumps in the tidal arms where 
the density has local maxima. An indication for the presence 
of such clumps in tidal arms can already be found in the ob- 
servations of Palomar 5 (Odenkirchen et al. 2001, 2003). The 
clumps have been first noticed in simulations by Capuzzo 
Dolcetta, di Matteo and Miocchi (2005) and were investi- 
gated further by di Matteo, Capuzzo Dolcetta & Miocchi 



Figure 8. A standing density wave has developed in the model 
K9. Top panel: At t = 0.87 Myr. A spherical cloud of tracer 
particles has been placed into the first clump (brown colored). 
Bottom panel: At t = 1.30 Myr (without the tracer particles). 

(2005). They noted already the wave- like nature of this phe- 
nomenon. 

The top panel of Figure [S] shows a spherical cloud of 
tracer particles (coloured brown). Figure [9] shows how the 
tracer particles have travelled further into the tidal arm 
while the position of the density maximum in the clump 
stayed (approximately) constant. Thus the clumps can be 
interpreted as wave knots of a density wave. 

A theoretical explanation for such clumps was published 
in Kiipper, Macleod & Heggie (2008). The top panel of Fig- 
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Figure 9. The model K9 at t = 1.7 A Myr. The tracer particles 
have travelled further while the wave maximum still persists. The 
leading and trailing arms have wound up. They are separated by 
the potential wall of the effective potential. The tip of the leading 
arm has hit the remnant of the star cluster again. We applied a 
weighting exponent to the particle mass in the expression (|B1||. 



ure|10|sho-ws for a few particles that they move on cycloids. 
The clumps appear at the position where many of the loops 
or turning points of the cycloids overlap. For a more detailed 
theory, see Just et al. (2009). The bottom panel of Figure [To| 
shows the radius as a function of time. We find approximate 
harmonic motion in the tidal arms. 

The dotted lines from the Galactic centre in Figures [7j 
[S] and [9] show the angles <po between density maxima in the 
leading and trailing arm. In order to plot these angles we 
determined the w coordinate of the maxima in the mean 
density (cf. the top panel of Figure 15 below) and obtained 



the corresponding Cartesian coordinates from our data files. 

Figure [Tl] shows the histogram of the epicentre radii 
Reo of the cycloid orbits, the epicyclic periods T^o, and 
the dimensionless angular momentum differences AL / Lc ~ 
(L — Lc)/Lc for different times, where Lc is the angu- 
lar momentum of the circular orbit. For the epicentre radii 
(and the epicyclic periods), stars within twice (and once) 
the membership radius were not included in the statistics. 
The epicentre radii are given by the arithmetic mean of the 
last maximum and minimum in the epicyclic amplitude. The 
epicyclic period is given by the time between the last two 
minima in the epicyclic amplitude. Note that ai t — 0.87 
Myr not all particles have completed one epicyclic period. 
For the dimensionless angular momentum differences we in- 
cluded only stars within 25 degrees around the density max- 
imum in the clump. Thus one can see two side lobes in all 
panels corresponding to the leading and trailing arms. Fig- 
ure |12| shows the distribution of Jacobi energies Ej^i and the 
histogram of peri- and apocentre radii of the cycloid orbits 
in the tidal arms. 



o 




0.0 0.5 1.0 1.5 2.0 2.5 

t [Myr] 



Figure 10. Top panel: Motion on cycloids within the tidal arms 
for a few particles in the model K9. The particles escape from 
locations near the cluster centre at (-20,0) either into the leading 
or the trailing arm. Bottom panel; Amplitude as a function of 
time for the same orbits. Small deviations from the harmonic 
motion can be seen which may be apparent deviations due to a 
slight change of the orbital frequency by dynamical friction or a 
deflection by 2-body encounters. 



The angle ifio can be expressed as 



10 A. Ernst, A. Just, R. Spurzem 





1.67 1.68 1.69 1.70 1.71 1.78 
Ej [10^ pcVMyr^] 



0.4 0.5 
T«„ [Myr] 




Leading 



Trailing 




0.4 0.6 



Figure 12. Top panel: Distribution of Jacobi energies. Bottom 
panel: Histogram of peri- and apocentre radii of the cycloid orbits 
in the tidal arms. 




Figure 11. Histograms for the leading and trailing arms. All par- 
ticles within twice (top panel) or once (middle panel) the initial 
tidal radius or outside of 25 degrees around the density maxi- 
mum in the first clump have been excluded from the statistics. 
We neglected deviations from the harmonic motion for the top 
and middle panel. Top panel: Epicentre radius distribution. Mid- 
dle panel: Epicyclic period distribution. We set the bin frequency 
equal to the data output frequency of our code in order to avoid 
an unphysical higher harmonic in the histogram. Bottom panel: 
Distribution of dimensionless angular momentum differences. 



Figure 13. The z-component of angular momentum with respect 
to the Galactic centre for a few orbits. In the cluster the angular 
momentum is not conserved since the cluster potential breaks 
the axisymmetry of the effective potential. In the tidal arms the 
angular momentum is only approximately conserved. 




where ujc and u) are the circular frequencies at the radius 
Rc of the circular orbit and in the vicinity of Rc , P is given 
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by Equation (271 and l\L/Lc is the most frequent dimen- 



sionless angular momentum difference. 

The subscript "L" refers in the following discussion to 
quantities which are expressed as a function of AL/Lc|^ 
while the subscript "0" refers to quantities which are directly 
measured from the simulation. 

The epicentre radius of the cycloids is given by 



Rel = i?c 1 + 



(37) 



According to Just et al. (2009), the epicyclic amplitude 
can be expressed as 



3-Q 



AL2 ASj 



2a + 2'"'^ Ll 



+ 



(38) 



to second order in the dimensionless angular momentum dif- 
ference, where A_Bj = — $off,tid(J?c) is the Jacobi energy 
difference with respect to the effective tidal potential at Rc ■ 
Table [3] shows a comparison of the measured angles, 
epicentre radii, epicyclic amplitudes, peri- and apocentres 
and the theoretical estimates from the dimensionless angular 
momentum difference. We give the measured angle </5o, the 
estimate ipL according to Equation (36 1, the error /^.tp/tpo — 
{ipo — fL)/'Po in percent, the most frequent epicentre ra- 
dius Reo in Figure [Tl] the epicentre radius Rel according 
to Equation (37 1, the most frequent epicyclic period T^o in 
Figure f 1 and the theoretical epicyclic period T^[Rel) using 



Equation (26 1 with the epicentre radius Rel- Furthermore, 
we give the tidal radius n = {GM^i/ P'^)ujI]Y^^ 
(King 1962), the arc length yo — Rc^po, the A factors (Just 
et al. 2009), 



A. 



yO 



1 Va + 1 yo 



a rt 



\Re,l-Rc\ 
rt 



(39) 



where AyQ is a first-order approximation and the error 
AA/Ayo ~ [Ayo — AL)/Ayo is in percent. We also give the 
most frequent peri- and apocentre radii Rpo and Rao from 
Figure [T2| the most fr equ ent scaled Jacobi energy difference 
^Ej juj^ from Figure 12 the epicyclic amplitude ArL from 
Equation ( 38 1 and the obtained peri- and apocentre radii 
RpL and Ral, where \Rpl - Rel\ = \Ral - Rel\ = Ar^,. 

There are systematic errors in both $l and Al and 
also in _RpL and Ral- The reason is shown in Figure [Ts] In 
the tidal arms the angular momentum is only approximately 
conserved. The reason is the influence of the cluster poten- 
tial which breaks the axisymmetry of the effective potential. 
However, in the cluster the angular momentum changes in 
a much shorter time scale. An estimate for the cumulative 
perturbation AL of L is given by 



AL = 



I Jin X a) dt| 



(40) 



Vdrift 

where a, R= [Rj, $ci and Vdrift are the acceleration, galac- 
tocentric radius, potential energy of the cluster and the drift 

* except for the case of the tidal radius r^, where the subscript 
"L" refers to the Lagrangian points Li and L2 
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Figure 14. Top panel: Time evolution of the angular momentum 
difference between the cluster center and the position in the tidal 
arm for the model K9. The angular momentum difference is nor- 
malized by the angular momentum of the cluster center. In order 
to show the asymmetry between the leading and trailing arms, 
the lines for the leading arm have been rotated by 180 degrees 
about the origin and replotted in grey. Bottom panel: As in the 
top panel, but for the energy (internal and external) difference 
between the cluster center and the position in the tidal arm for 
the model K9. The energy difference is normalized by the energy 
of the cluster center. 



velocity, respectively. Thus a slow drift velocity increases the 
change in L. Here a more detailed theory is desirable. 

For Figures [l4] - 116[ we averaged over spheres with a 
radius which was approximately equal to the width of the 
tidal arms in the xy plane (see Appendix [b] for the details). 
Since the Figures were still noisy, we have used a median 
smoothing in addition. The width of the smoothing kernel 
has been taken to be twice the membership radius defined 
in Equation (|3|. 

For four different times, the top panel of Figure [14] 
shows the z-component of the dimensionless mean angular 
momentum difference [Iziw) — lz{0)] /lz{0) along the tidal 
arms. The specific angular momentum was calculated with 
respect to the Galactic center. In order to show the asym- 
metry between the leading and trailing arms, the lines for 
the leading arm have been rotated by 180 degrees about the 
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Table 3. Comparison of measurement and theory for the angles of the density maxima and the A factors. For explanations see the text. 



origin and replotted in grey. This kind of asymmetries arise 
from the geometry of the effective potential. 

The bottom panel of Figure [l4] shows the same for the 
mean energy difference [e{w) — e(0)] /e(0). The specific en- 
ergy was calculated with respect to the Galactic center. For 
the definition of the energy, see Section [2.3| A positive en- 
ergy difference corresponds to the trailing arm while a neg- 
ative energy difference corresponds to the leading arm. This 
is in accordance with the positive normalization in Equation 
( |25[ ) for the scale free potential in Equation \21\ . This Figure 
also shows an asymmetry between leading and trailing arm. 
What can be seen in this plot is that with the time more 
and more particles with a low energy difference with respect 
to the cluster center stream into the tidal arms. Thus the 
modulus of the mean energy difference falls off with time. 
In this connection it is worthwhile to mention that the par- 
ticles with a low Jacobi energy stream into the tips of the 
tidal arms. This can be seen in Figure |9] The stars in the 
tips of the tidal arms are far away from the potential wall 
of the effective tidal potential which lies below the solid line 
of the orbit of the star cluster center. 

The top panel of Figure [15] shows the density profile 
along the tidal arm coordinate w. At t = 0.43 Myr, one 
clump can be seen in the leading arm. At t = 0.87 Myr, 
two clumps can be seen in the leading arm and one in the 



trailing arm. At t = 1.30 Myr, three clumps can be seen in 
the leading arm and two in the trailing arm. 

The bottom panel of Figure [15] shows the profile of the 
ID velocity dispersion along the tidal arm coordinate w. The 
averaging for the calculation of the velocity dispersion has 
been done in small spheres around the tidal arm coordinate 
system whose radius was approximately 1/5 of the width of 
the tidal arms in the xy plane (see Appendix [b] for more 
details). The velocity dispersion profile also exhibits local 
maxima at the positions of the density maxima. This is in 
accordance with the notion that the clumps in the tidal arms 
occur at the positions where many of the loops or turning 
points of the cycloid orbits overlap. At these positions, the 
random velocities should exhibit maxima as well. Note that 
the first maximum in the leading arm at t = 0.43 Myr cannot 
yet be seen clearly in the bottom panel of Figure [15] Figure 
[7] shows that this density maximum is still in the process of 
building up. 

Figure [16] shows the profile of the cluster gravitational 
potential along the tidal arm coordinate w at four different 
times. One can see that the potential well of the star cluster 
is deeper in the beginning but gets shallower as the cluster 
loses mass. 

Figure [it] shows the evolution of the cluster mass con- 
tained within the tidal radius for the model K9. In addition, 
the mass in the tidal arms is shown as a function of time. 
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Figure 15. Top panel: Time evolution of the mean density of 
stars along the tidal arms for the model K9. One can see that 
several density wave maxima develop with time. Bottom panel: 
Time evolution of the ID velocity dispersion along the tidal arms 
for the model K9. The characteristic features in the mean density 
(top panel) can also be seen in the velocity dispersion. 
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Figure 17. Evolution of the cluster mass within the tidal radius 
and the mass in the tidal arms for the model K9. It can be seen 
that more particles escape into the trailing arm than into the 
leading arm. The ratio of leading arm mass to trailing arm mass 
is always roughly 85%. The thin vertical dotted lines correspond 
to t = 0.22, 0.43, 0.87, 1.30 and 1.74 Myr. 



It can be seen that more particles escape into the trailing 
arm than into the leading arm. In the relaxation-driven dis- 
solution scenario this would be paradoxical since the inner 
Lagrange point Li is at a lower energy than the outer La- 
grange point L2 according to Figure [4] However, most stars 
are in the high-energy region of the star cluster. For these 
particles the phase space for escape into the trailing arm is 
larger than that for escape into the leading arm. 



Leading arm 
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Figure 16. Time evolution of the mean star cluster potential 
along the tidal arms for the model K9. One can see that the 
potential well is deeper in the beginning but gets shallower as the 
cluster loses mass. 



4.2 Lifetime scaling and RE classification 

Figure [I8] shows the scaling of the half-mass time as a func- 
tion of the particle number. The half-mass time is the time 
after which the star cluster has lost half of its initial mass 
due to escaping stars. It also shows as dots the times, when 
the cluster has one third or one fifth, respectively, of its 
initial mass. The particle number in Figure [18] ranges from 
iV = 10^ npto N = IQ^. The fit of a power law shows 
that the half-mass time depends only slightly on the particle 
number A'^. 

This is in contrast to the relaxation-driven dissolution 
of star clusters. If the dissolution is relaxation-driven, the 
stars are scattered above the escape energy (or the critical 
Jacobi energy) by two-body relaxation before they can es- 
cape through exits in the equipotential surfaces around the 
Lagrangian points Li and L2- In this case the half-mass or 
dissolution time should depend more strongly on the parti- 



cle number than in Figure 18 Baumgardt (2001) developed 
a detailed theory for the relaxation-driven dissolution of star 
clusters on circular orbits in a steady tidal field with back- 
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Figure 18. Scaling of the ratio of half-mass time ty^^l f to crossing 
time tcr as a function of the particle number TV for the models 
Kl - K9. The errors correspond to 1 crossing time tcr- We have 
tcr = 0.123 Myr in all models. The "third-" and "fifth-mass" 
times are also shown. 
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Figure 19. Projection of initial Jacobi energies of stars onto the 
a;— "iE"eff plane of Figure[4]for the model K9. The dashed lines mark 
the different regions of the Radius-Energy (RE) classification. 



scattering of potential escapers in which the half-mass time 
scales as thaif oc t%*. 

Fukushige & Heggie (2000) give a hint for the under- 
standing of the scaling of the half-mass time in our mod- 
els. In Figure [19] one can see that the cluster fills the en- 
ergetic region above the total effective potential of Fig- 
ure |4] Many stars are initially outside of the tidal radius 
tl ~ [x{Li) + x{L2)] /2 and most particles have Jacobi en- 
ergies Ej per unit mass which are higher than the mean 
effective potential Ej,l = [Ej{Li) + Ej{L2)] /2 of both La- 



grange points Li and L2. We initially have for the model K9 
the following ratio of particle numbers: 



§^^0.37 and ^^^^ 



(41) 



The stars which are outside of tl and the high-energy parti- 
cles with respect to the critical Jacobi energy Ej^l can leave 
the cluster relatively fast as compared with the relaxation 
time, provided they are not bound by a non-classical inte- 
gral of motion which would hinder their escape. Figure [S] 
shows that above a certain Jacobi energy threshold in the 
high-energy regions all orbits are chaotic and not subject to 
a non-classical integral of motion. 

It is possible to classify the particles initially according 
to their membership to one of four regions: 

(i) Large Radius High Energy (LRHE) region 

(ii) Small Radius High Energy (SRHE) region 

(iii) Large Radius Low Energy (LRLE) region 

(iv) Small Radius Low Energy (SRLE) region 

The distinction between these regions is shown with dashed 



lines in Figure 19 We call this the Radius-Energy (RE) clas- 



sification. The classification arises due to the existence of the 
Lagrange points Li and L2 at proximate (or equal) Jacobi 
energies. For the model K9 with — 10^, we initially have 
the following occupation numbers of the four regions, 



iVLRHE = 24370, iVsRHE = 39854, 
iVLRLE = 2831, A^sRLE = 32945. 



(42) 



In an exact treatment the critical equipotential line should 
be taken as the dividing line between small and large radii in 
the RE classification. However, in real Ai'-body simulations 
it is more convenient to adopt the tidal radius for this pur- 



pose as has been done in the counting for Equations ( 42 1 



Particles in the LRLE region are immediately lost due to 
the energy barrier if they are not bound to the cluster by a 
non-classical integral of motion. Particles in the LRHE and 
SRHE regions can mix dependent on their individual posi- 
tion and velocity. Particles in the SRLE region are bound 
to the cluster until they are lifted to the SRHE region by 
secular evolution or 2-body relaxation. 
Particularly the ratio 



(43) 



A^LRLE + .M'lRHE + AfsRHE 
AfsRLE 

determines the physics of the dissolution process, where 
Mlrle, Mlrhe, Msrhe and Msrle are the occupation 
masses of the four regions. If om is close to zero the main 
process which leads to the dissolution of the cluster is two- 
body relaxation, which scatters stars from the SRLE region 
into the two high-energy regions. The larger the particle 
number A'' is, the slower is this process. We speculate that 
qa/ was very small for the old globular clusters in the halo 
of the Milky Way and that their dissolution is relaxation- 
driven, but that many young star clusters (open clusters 
or associations) with larger values of oa/ may form at all 
times in the Milky Way and dissolve fast as compared with 
the Hubble time. If om is sufficiently large, mass loss from 
the SRLE region seems to be dominated by a self-regulating 
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process of increasing Jacobi energy due to the weakening of 
the potential well of the star cluster, which is induced by 
the mass loss itself (Just et al. 2009). A simple estimation 
shows that the critical Jacobi energy Ej^l increases more 
slowly with time as compared with the Jacobi energy Ej of 
a star in the non-stationary gravitational potential of the 
star cluster. While the particles in the LRLE, LRHE and 
SRHE regions of the star cluster move away from the clus- 
ter along the tidal arms, particles are continually shifted 
from the SRLE region into the two high-energy regions as 
the potential well of the star cluster gets shallower (cf. Fig- 
ure 161. In addition, a fraction of particles is scattered from 



the SRLE region into the high-energy regions by two-body 
relaxation. The two-body relaxation leads to the small slope 
0.127 ±0.030 in Figure It is small since qa/ is very large 
for the models Kl - K9. From the values in Equation (42 1 we 
obtain a^/ ~ 2 for the model K9 with the valid assumption 
that the particles of different mass are initially uniformly 
mixed in radius and Jacobi energy per unit mass. 

If the physical tidal radius is equal to the radius where 
the density of the star cluster (King) model vanishes, we 
have A'^LRLE = Alrhe = and only two of the four regions 
are occupied with particles. This is the standard case used 
in A-body simulations of star clusters in a tidal field so 
far (e.g. Baumgardt & Makino 2003, Trenti, Heggie & Hut 
2007, Ernst et al. 2007). In this case, qm is small (typically 
a few percent) and the dissolution times are considerably A- 
dependent as we checked with a few models (A — 10^,2 x 
10^5 X 10^ 10^,3 X 10*) using nbody6gc. Furthermore, in 
this case our preliminary models suggest that the dissolution 
time directly scales with a power of the relaxation time. 
However, a more detailed study seems to be essential. On the 
other hand, Tanikawa & Fukushige (2005) adopted initial 
models where the King cutoff radius was smaller or larger 
than the physical tidal radius. By decreasing the size of the 
initial star cluster further am can be forced to vanish. 

We argue that the situation that the cluster is divided 
into the four regions of the RE classification (with certain 
occupation numbers and masses) is the typical situation 
for newly formed star clusters. A first crucial question is 
whether stars can form in all regions. The answer is yes, if 
the condition for star formation is fulfilled. According to the 
modern picture of gravo-turbulent star formation (e.g. Mac 
Low & Klessen 2004, Ballesteros-Paredes et al. 2007), super- 
sonic turbulence and shocks create initial density enhance- 
ments in a molecular cloud. The formed molecular cloud core 
contracts gravitationally and fragments eventually. Finally, 
protostellar seeds form, accrete in-falling material and be- 
come main sequence stars. The condition for star formation 
is independent of the distinction between high- and low- 
energy regions of the effective potential. Thus one would ex- 
pect that stars form initially in the high-energy regions and 
the SRLE region slowly builds up as more material moves 
towards the center of the new star cluster. Due to the turbu- 
lent structure within the molecular cloud it is also possible 
that a small fraction of stars forms in the LRLE region. 

In the Galactic center, the supersonic shock and turbu- 
lent velocities must be high enough to form mean densities 
which withstand the tidal shear forces. According to Morris 
(1993), the critical mean number density for gravitationally 
bound clouds in the Galactic center region is given by 
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ricrit = 10 cm 



1.6 pc 
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(44) 



where Rg is the galactocentric radius. 

The picture sketched above would be similar if the star 
cluster formation in the Galactic center is triggered by the 
collision of two clouds. For typical parameters (e.g. for the 
formation of clusters like Arches and Quintuplet) the rate 
of such cloud collisions in the Galactic center is low as com- 
pared with the reciprocal of the lifetime of OB stars and can 
be crudely estimated to be 
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where Mdoud, and are the mass, the column density 
and the velocity dispersion of a cloud (Hasegawa et al. 1994, 
Stolte et al. 2008). 

Finally, we note that the Jeans time scale is of the same 
order as the dissolution times of our models in the Galac- 
tic center. According to Hartmann (2002), who explored an 
earlier idea by Larson (1985), the Jeans (or fragmentation) 
time scale of a gaseous filament can be written as 
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10 K 



1/2 



Myr 



(46) 



where T is the temperature and Ay « 5 is the visual extinc- 
tion through the center of the filament (see also Klessen et 
al. 2004). 



5 DISCUSSION 

We have studied the dissolution of star clusters in an ana- 
lytic background potential of the Galactic centre by means 
of direct A-body simulations. We described in detail the al- 
gorithm of our parallel A-body program nbodyOgc which 
is based on Aarseth's series of A-body codes (Aarseth 1999, 
2003, Spurzem 1999). It includes a realistic dynamical fric- 
tion force with a variable Coulomb logarithm based on the 
paper by Just & Penarrubia (2005) . The initial value for the 
circular orbit of the model K9 is In A ~ 1.7. It turns out 
that, even for a 10® Mq cluster, the dynamical friction force 
is too weak to let a cluster on a circular orbit at Rc = 20 
pc spiral into the Galactic centre within the lifetime of its 
most massive stars. Thus we did not resolve the "paradox 
of youth" (Ghez et al. 2003). 

However, we have studied in detail the dynamics of dis- 
solving star clusters on circular orbits in the Galactic centre. 
The key to the understanding of this dynamical problem is 
the gravitational potential which is the sum of the effective 
tidal potential and the star cluster potential. Along the orbit 
of the star cluster, the effective tidal potential resembles a 
parabolic wall. However, in the close vicinity of the Galactic 
centre there are deviations from the parabolic shape due to 
higher-order terms in the Taylor expansion of the effective 
tidal potential. Due to this asymmetry, the Lagrange points 
Li and L2 lie at different energies. 
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We have studied in detail the properties of the tidal 
arms of a dissolving star cluster in a galactic centre. The 
density wave phenomenon found by Capuzzo Dolzetta, di 
Matteo & Miocchi (2005) and di Matteo, Capuzzo Dolcetta 
& Miocchi (2005) appears in our model K9. The angles of 
the clumps can be calculated with the theory from Just et 
al. (2009). 

We have presented a method to study the structure 
of tidal arms by using an eigensolver. The eigensolver 
calculates numerically a ID coordinate system along the 
tidal arms and calculates characteristic dynamical quanti- 
ties along this coordinate system. 

It may be of interest to note that more particles escape 
into the trailing tidal arm than into the leading tidal arm. 
For the high-energy particles the phase space for escape into 
the trailing arm is larger than that for escape into the leading 
arm. The fractions of initial conditions in the phase space 
for escape into the leading and trailing arm, respectively, 
have to be calculated numerically for several Jacobi energies. 
The result would be called the 'basins of escape' of the star 
cluster (e.g. Aguirre et al. 2001, Contopoulos 2002, Ernst et 
al. 2008). This kind of asymmetry between the arms does 
not depend on the particle number. 

The half-mass times of our models Kl - K9 depend only 
weakly on the particle number which indicates that two- 
body relaxation is not the dominant mechanism leading to 
the dissolution. The reason is that the initial models are di- 
vided into four different regions in radius and specific Jacobi 
energy space. This division has been termed the Radius- 
Energy (RE) classification. The division of a newly formed 
star cluster into the four regions of the RE classification 
is probably a typical situation which is consistent with the 
modern picture of gravoturbulent star formation (e.g. Mac 
Low & Klessen 2004, Ballesteros-Paredes et al. 2007). If the 
ratio um (which has been defined in Section 4.21 is large 
enough, the dissolution is no longer relaxation-driven but 
the mass loss from the SRLE region is governed by a self- 
regulating process of increasing Jacobi energy due to the 
weakening of the potential well of the star cluster, which is 
induced by the mass loss itself (Just et al. 2009). Predictions 
about the fractions of stars which belong to the four different 
regions (i.e., the occupation numbers and masses) may be 
an important result of the emerging theory of star cluster 
formation. What are typical ratios of occupation numbers 
and masses in regions with efficient star formation? How 
do the occupation numbers, masses and their ratios differ 
between open and globular clusters? From the side of stel- 
lar dynamics the scaling problem of the dissolution times 
(Baumgardt 2001) needs to be solved for the new dissolution 
mechanism due to a non-stationary gravitational potential 
combined with the effect of two-body relaxation. In this pa- 
per, we also conjecture that the ratio qm can be used to 
draw a distinction between associations, open and globular 
clusters, i.e. that the old globular clusters obey q:a/ ^ 1 and 
that their dissolution is relaxation-driven, while the open 
clusters and associations obey aM ^ 0. 
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Figure Al. Taylor expansion of the effective tidal potential for 
the scale free model up to the 5th order. 



APPENDIX A: TAYLOR EXPANSION OF THE 
EFFECTIVE TIDAL POTENTIAL 

Usually a Cartesian Taylor expansion of the effective tidal 
potential is used. Up to 5th order, the 3D Cartesian Taylor 
expansion of the effective tidal potential for the scale free 
model is given by 
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(Al) 



where _Ro and uoq are the radius and the frequency of the 



circular orbit. This solution is shown in Figure Al (which 
may be compared with Figure [2|. It can be seen that this 
Taylor expansion cannot be used to study the properties of 
tidal arms in the Galactic centre. For extended tidal tails 
cylindrical coordinates should be used and for the radial 
asymmetry the exact potential. 

The expansion up to the second order coincides with 
the tidal approximation. We have for the scale free model 
(a — 3)a;o/2 = (kq ~ 4tJo)/2 which is the coefficient of the 
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second-order term in the tidal approximation, where kq is 
the epicychc frequency. 



APPENDIX B: 
SYSTEM 



TIDAL ARM COORDINATE 



Based on routines from Numerical recipes (NR, Press et 
al. 2001), we developed the eigensolver eigentid which 
calculates numerically a ID coordinate system along the 
tidal arms and evaluates characteristic dynamical quanti- 
ties along this coordinate system. We denote the coordinate 
along the tidal arms as w, where negative values refer to 
the leading arm and positive values to the trailing arm. The 
NR routine tred2 uses the Householder reduction of a real 
symmetric n x n matrix to convert it to a tridiagonal form. 
The NR routine tqli uses the QL algorithm to determine 
the eigenvalues and eigenvectors of the matrix which has 
been brought into tridiagonal form before (see NR, Chap- 
ters 11.2 and 11.3). We use the tensor of inertia and denote 
the eigenvectors corresponding to the minimum eigenvalue 
c, the medium eigenvalue b and the maximum eigenvalue a 
as the minimum, medium and maximum eigenvectors, re- 
spectively. Then the algorithm proceeds as follows: 

(i) Read snapshot with particle masses, positions and ve- 
locities in the cluster rest frame. 

(ii) Calculate optionally gravitational potential and den- 
sity (using the method by Casertano & Hut 1985) for this 
snapshot. 

(iii) Start calculation in the origin of coordinates (0, 0, 0). 

(iv) Obtain a neighbour sphere with radius i?cut and cal- 
culate its centre of mass {xcm,ycm, Zcm)- 

(v) Calculate physical quantities averaged over the neigh- 
bour sphere: Mean specific angular momentum, mean spe- 
cific energy, mean density, velocity dispersion, mean poten- 
tial. Write all quantities to a data file. 

(vi) Calculate the tensor of inertia of the neighbour 
sphere with respect to the centre of mass of the neighbour 
sphere. It is given by 



We remark that the inertia ellipsoids of the neighbour 
spheres have an oblate shape, i.e. the three eigenvalues a, b, c 
of Qjk satisfy a ~ b > c. Also, a weighting exponent can be 
assigned to the particle mass in the expression (Bl I. In this 



case, the eigensolver follows the mass distribution within the 
tidal arms in a different way. This method has been applied 
for Figure [9] 

This paper has been typeset from a TjnjX/ M^rjX file prepared 
by the author. 
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where Axi = Xi — x^m , Ayi ^ yi — ycm and Azi = Zi — z^m 
are the relative positions of the ith particle in the neighbour 
sphere with respect to its centre of mass and rrii is the mass 
of the ith particle. 

(vii) Calculate the eigenvalues and eigenvectors of Qjk- 

(viii) Go along the direction of the maximum/medium 
eigenvector until a critical density is reached to find the new 

Rcut- 

(ix) Check for acute angle between previous and current 
minimum eigenvector. If the angle is acute, change the sign 
of the eigenvector. 

(x) Go one step along the direction of the minimum eigen- 
vector. 

(xi) Repeat from 4. until the particle number within the 
neighbour sphere drops below a certain threshold as the first 
tidal arm ends. 

(xii) Start from 3. for the second tidal arm. 



